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Abstract. We consider the symmetric Darlington synthesis of a p x p rational symmetric Schur 
function S with the constraint that the extension is of size 2p x 2p. Under the assumption that S is 
strictly contractive in at least one point of the imaginary axis, we determine the minimal McMillan 
degree of the extension. In particular, we show that it is generically given by the number of zeros 
of odd multiplicity of I p — SS* . A constructive characterization of all such extensions is provided 
in terms of a symmetric realization of S and of the outer spectral factor of I p — SS*. The authors's 
motivation for the problem stems from Surface Acoustic Wave filters where physical constraints on 
the electro-acoustic scattering matrix naturally raise this mathematical issue. 
Keywords, symmetric Darlington synthesis, inner extension, MacMillan degree, Riccati equation, 
symmetric Potapov factorization. 

1 Introduction 

The Darlington synthesis problem has a long history which goes back to the time when computers 
were not available and the synthesis of non-lossless circuits was a hard problem: the brilliance of the 
Darligton synthesis was that it reduced any synthesis problem to a lossless one. In mathematical 
terms, given a(pxp) Schur function S, say, in the right half-plane, the problem is to imbed S into 
a (m + p) x (m + p)-inner function S so that: 

5 = ( Sn S <S ) ' S(iu)S*(icj) = I m+P , wGR. (1) 

This problem was first studied by Darlington in the case of a scalar rational S [Dlj , was generalized 
to the matrix case [Bl] , and finally carried over to non-rational S \A2\ ID2j . We refer the reader 
to the nice surveys [BM} fD3] for further references and generalizations (e.g. to the non-stationary 
case). An imbedding of the form (P) will be called a Darlington synthesis or inner extension, or 
even sometimes a lossless extension of S. 

A Darlington synthesis exists provided that S(iuj) has constant rank a.e., that the determinant of 
I p — S(iu)S* (ioo) (viewed as an operator on its range) satisfies the Szego condition (see e.g. [G]). 
and that S is pseudo-continuable across the imaginary axis meaning that there is a meromorphic 
function in the left half-plane whose nontangential limits on iM. agree with S(iu) a.e. [A2} IA31 ID2} 
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D3, DHJ. If moreover 5 has the conjugate-symmetry, then S may be chosen with this property. 
When S is rational the above conditions are fulfilled so that a Darlington synthesis always exists. 
In addition, it can be carried-out without augmenting the McMillan degree (i.e. we may require in 
(PQ) that deg5 = degS) and choosing m = p; this follows easily from Fuhrmann's realization theory 
[Fl| and the arguments in [D3j or else from more direct computations carried out in \AV\ IGR| . In 
particular S can be chosen rational, and also to have real coefficients if S does. 
When S is the scattering matrix of an electric p-pole without gyrators |Bl], the reciprocity law 
entails that S is symmetric and the question arises whether the extension S can also be made 
symmetric; this would result in a Darlington synthesis which is itself free from gyrators. In [AV^ 
it is shown that a symmetric Darlington synthesis of a symmetric rational S indeed exists and, 
although one can no longer preserve the degree while keeping m = p, he can at least ensure that 
deg5 < 2degS. The existence of a symmetric Darlington synthesis for non-rational functions has 
been studied in |A3j, in the slightly different but equivalent setting of J-inner extensions. 
In |AV| it is also shown that, by increasing the size m to p + n, where n is the degree of S, it 
is possible to construct a symmetric extension of exact degree n. However, such an increase of m 
is not always appropriate. In fact, although the original motivations from circuit synthesis that 
brought the problem of lossless imbedding to the fore are mostly forgotten today, the authors of the 
present paper were led to raise the above issue in connection with the modeling of Surface Acoustic 
Waves filters [BEGO]. In this context, physical constraints impose m = p, so that each block of 
the electro-acoustic scattering matrix S in ([1]) has to be of size p x p. 

It is thus natural to ask the following : given a symmetric rational S, what is the minimal degree 
of a symmetric lossless extension S? This is the problem that we consider. For scalar systems, 
this minimal degree has been known for a while and can be found, for instance in |YWPj (see also 
Section 3 below). The present paper will generalize this to the matrix- valued case. We restrict our 
attention to the case where S is strictly contractive in at least one point of the imaginary axis. 
This implies that the extension will have size 2p. For the general case, that is, with extensions of 
lower size, the analysis seems to be more difficult and it will possibly be treated in a subsequent 
paper. 

In Section [2] we introduce some notations. In Section [31 we shed light on the problem by discussing 
the elementary scalar case, that is, p = 1. In Section [4] we recall some results of Gohberg and 
Rubinstein [GR] about a state space construction of an inner extension preserving the degree and 
we characterize all inner extensions in terms of minimal ones. In Section [5] we present a simple 
method to construct (possibly unstable) symmetric unitary extensions. In Section [6] we finally 
produce a symmetric inner extension of minimal degree. In Section [7] we discuss the symmetric and 
conjugate symmetric unitary extension of a rational symmetric Schur function which is conjugate 
symmetric (i.e. that has real coefficients), and we show on an example that its minimal degree is 
generally larger than the one attainable without the conjugate-symmetry requirement. 

2 Preliminaries and notations 

Throughout, if M is a complex matrix, we let Tr(M) stand for its trace, M T for its transpose, and 
M* for its transpose-conjugate. We denote respectively by 

n+ = {s £ C; Re s > 0} and IT = {s e C; Re s < 0} 

the right and left half-planes. 
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In System Theory, a rational function whose poles lie in n~ is called stable, and a rational function 
which is finite (resp. vanishing) at infinity is called proper (resp. strictly proper). System Theory 
is often concerned with functions having the conjugate symmetry: W(s) = W(s), but we shall not 
make this restriction unless otherwise stated. A rational function has the conjugate-symmetry if, 
and only if, it has real coefficients. 

For W(s) a matrix-valued function on C, we define its para-hermitian conjugate W* to be: 

W*{s) := W(-s)*. (2) 

Note that * has two different meanings depending on its position with respect to the variable; this 

slight ambiguity is common in the literature and allows for a simpler notation. 

Note that W*(iuj) = W(iuj)* on the imaginary axis, and if W is a polynomial then W* is also a 

polynomial whose zeros are reflected from those of W across the imaginary axis. 

We say that a rational p x m matrix-valued function S holomorphic on II + is a Schur function if 

it is contractive: 

S(s)S(s)* <I P , s G n+. (3) 
A rational p x p Schur function S is said to be lossless, or inner, if 

S(wj)S(icj)* = I P , uj G M. (4) 

The scalar rational inner functions are of the form q* jq where q is a polynomial whose roots lie in 
IT~; if degg = n, such a function is called a Blaschke product of degree n. A (normalized) Blaschke 
product of degree 1 is just a Mobius transform of the type 

h(s):=(s-0/(s + 0, £GII+. (5) 

The natural extension to the matrix case is given by 

which is the most elementary example of an inner function of degree 1. Actually, it is a result of 
Potapov [P |ID4| that these and unitary matrices together generate all rational inner matrices. More 
precisely, if Q is such a matrix and £i, . . . its zeros (i.e. the zeros of its determinant which is a 
Blaschke product) ordered arbitrarily counting multiplicities, there exist complex unitary matrices 
Ui, . . . , U n+ \ such that 

Q = U 1 % U 2 % . . . U n % U n+1 . (7) 

An inner matrix like U\ U2 is often called an elementary Blaschke factor. 
Given a proper rational matrix S, we shall write 



S 



A 


B 


C 


D 



(8) 



whenever (A, B,C, D) is a realization of S, in other words whenever S(s) = C(sl n — A)~ 1 B + D 
where A, B, C, D are complex matrices of appropriate sizes. Because S in this case is the so-called 
transfer function [AV] IBGKl IKFA| of the linear dynamical system: 

x = Ax + Bu, y = Cx + Du, (9) 
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with state x, input u, and output y, we say sometimes that A is a dynamics matrix for S. The 
matrices C and B are respectively called the output and input matrices of the realization. 
Every proper rational matrix has infinitely many realizations, and a realization is called minimal 
if A has minimal size. This minimal size will be taken as definition of the McMillan degree of S, 
abbreviated as degS. As is well-known [KFA, BGK], the realization ([8]) is minimal if and only if 
Kalman's criterion is satisfied, that is if the two matrices: 

[B AB ... A n ~ 1 B } , [C T A T C T . . . (A T ) n ~ l C T ] , (10) 

are surjective, where n denotes the size of A. 

The surjectivity of the first matrix expresses the reachability of the system, and that of the second 
matrix its observability. Any two minimal realizations can be deduced from each other by a linear 
change of coordinates: 

(A, B, C, D) ^ (TAT^^B, CT~ l ,D), T an invertible matrix, 

so that a dynamics matrix of minimal size for S is well-defined up to similarity. In particular the 
eigenvalues of A depend only on S and they are in fact its poles, the multiplicity of a pole being 
its total multiplicity as an eigenvalue by definition. The sizes of the Jordan blocks associated to 
an eigenvalue are called the partial multiplicities of that eigenvalue. The partial multiplicities may 
be computed as follows. Performing elementary row and column operations on S, one can put it 
in local Smith form at £ (see e.g. fGLRl sec.7.2.] or [BQ21IBGR] ): 

S(s) = £( S )diag[( S - V1 , (a " • ••,(«- 0" k , 0, • • • , 0]F(s) (11) 

where E(s) and F(s) are rational matrix functions that are finite and invertible at £ while k is the 
rank of S as a rational matrix and v\ < v 2 < . . . < are relative integers. These integers are 
uniquely determined by S and sometimes called its partial multiplicities at £. Note that £ is a pole 
if, and only if, there is at least one negative partial multiplicity at £. In fact, the negative partial 
multiplicities at £ are precisely the partial multiplicities of £ as a pole of S. 

One says that £ is a zero of S if the local Smith forms exhibits at least one positive partial 
multiplicity at £, and the positive partial multiplicities at £ are by definition the partial multiplicities 
of £ as a zero. If S is invertible as a rational matrix, it is clear from (jlip that the poles of S^ 1 are 
the zeros of S, with corresponding partial multiplicities. Note also that a zero may well be at the 
same time a pole, which causes many of the difficulties in the analysis of matrix valued functions. 
When 5 is inner, which is our main concern here, this does not happen because its poles lie in II - 
and its zeros in Il + . 

A rational matrix has real coefficients if, and only if, there exists a minimal realization which is 
real, i.e. such that A, B, C, and D are real matrices. As is customary in System Theory, we 
occasionally refer to a proper rational matrix as being a transfer function. If it happens to have 
the conjugate-symmetry, we say it is a real transfer function. 

The system-theoretic interpretation ([9]) of (|8|) makes it easy to compute a realization for a product 
of transfer- functions. In fact, if S\ is m x k and S2 is k x p, and if 



Si 



At 


B 1 ' 


Ci 





s 2 



A 2 


B 2 


c 2 


D 2 



(12) 
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then a short computation shows that the following two realizations hold: 



S 2 Si 














B 2 D 1 


D 2 d 


c 2 


D 2 D X t 



S 2 Si 



A 2 


B 2 C\ 


B 2 D X 







Bt 


c 2 


D 2 d 


D 2 Di j 



(13) 



Likewise, if k 
then 



m and D\ is invertible (so that S\ is a fortiori invertible as a rational matrix), 





f A 1 - B X D^C X 


B\D^ 1 






DI 1 



(14) 



Using (jlOp . it is immediate that the minimality of (|12p implies that of (|14p . In contrast, the 
realizations (|13p need not be minimal even if the realizations (|12p are: pole-zero cancellations may 
occur in the product S 2 Si to the effect that the multiplicity of a pole may not be the sum of its 
multiplicity as a pole of £2 (i-e. an eigenvalue of A 2 ) and as pole of Si (i.e. an eigenvalues of A{). 
One instance where (|13j) is minimal occurs when Si, S 2 have full rank and no zero of Si is a pole 
of S 2 and no zero of S 2 is a pole of Si (see [CN| ). This, in particular, is satisfied when Si and 
S 2 are inner, implying that the McMillan degree of a product of inner functions is the sum of the 
McMillan degrees. By ([7]), this in turn implies that the McMillan degree of an inner function is 
also the degree of its determinant viewed as a scalar Blaschke product. 

Whenever S is a transfer function, its transpose S T clearly has the same McMillan degree as S. A 
square transfer function S is called symmetric if S = S T , and then a realization is called symmetric 
if A = A T , B T = C and D = D T . It is not too difficult to see that a transfer function is symmetric 
if, and only if, it has a minimal realization which is symmetric [FHJ. The latter may be complex 
even if S is a real transfer function. 



3 The case of a scalar Schur function. 

For getting an idea of the solution to our problem, we first consider the symmetric inner extension 
of a scalar rational Schur function to a 2 x 2 inner rational function, that is we assume momentarily 
that p = m = 1. This case has been considered in [YWP]. Put 

q 

where pi and q are coprime polynomials such that degjpi} < deg{q}, pi is not identically zero, 
|pi(i(j)| < |<?(iw)| for cu E R, and q has roots in the open left half-plane only. The McMillan degree 
of S is just the degree of q in this case. As the orthogonal space to a nonzero vector v = (a, b) T G C 2 
is spanned by (— b, a) T , it is easily checked that every rational inner extension S of S, when all its 
entries are written over a common denominator, say, dq where d is a stable monic polynomial, is 
of the form 

e i8 2 " 

1 

where 9%, 9 2 S R and p 2 is a polynomial solution of degree at most deg{dq} to the spectral factor- 
ization problem: 

dpi (dpi)* +p 2 p 2 = dq(dq)*, (15) 



dq 



1 



(dpi)* 

P2 



-Pi 

dpi 
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whose solvability is ensured by the contractivity of S. Clearly the extension is symmetric if and 
only if — which is compatible with (I15D if, and only if, all zeros of the polynomial 

dq(dq)* - dp\(dpi)* = dd*(qq* - pip\) 

have even multiplicity. Consider the polynomial 

M :=<?<?*- PlPi, ( 16 ) 

and single out its roots of even multiplicity by writing fi = (rir*) 2 r 2 r 2 , where r\ and r 2 are stable 
coprime polynomials and all the roots of r 2 are simple. For dd* \x to have roots of even multiplicity 
only, it is then necessary that r 2 divides d. Therefore, as the McMillan degree of an inner function 
is the degree of its determinant, we get 

deg{5} = deg{dq} > deg{r 2 <?}. (17) 

On another hand, a symmetric inner extension of McMillan degree deg{r2g} is explicitly given by 



q r 2 



ri r,r 2 



EL 



_£i 
q 



Pi 

q 



L2. o 

T2 

1 



(18) 



Thus we see already in the scalar case that the minimal attainable degree for a symmetric inner 
extension of S is the degree of S augmented by half the number of zeros of /x of odd multiplicity. 
Formulas (|17p and (|18p should be compared with the corresponding formulas (98) and (95) in 
[YWPj . 

As 1 — SS* = fJ-/(qq*), the zeros of \i are the zeros of 1 — SS* augmented by the common zeros to p\ 
and q* and the common zeros to p\ and q; the latter of course are reflected from the former across 
the imaginary axis, counting multiplicities. In particular a degree-preserving symmetric Darlington 
synthesis requires special conditions that can be rephrased as: 

(i) the zeros of 1 — SS* have even multiplicities, 

(ii) each common zero to S and (<S*) -1 , if any, is common with even multiplicity. 

Remark 3.1 Note that (i) is automatically fulfilled for those zeros located on the imaginary axis, 
if any, so the condition really bears on the non-purely imaginary zeros. Note also that (ii) concerns 
those zeros of S, if any, whose reflection across the imaginary axis is a pole of S; by the coprimeness 
of pi and q, such zeros are never purely imaginary. 



Our goal is to generalize the previous result to matrix-valued contractive rational functions. 



4 Inner extensions. 

We shall restrict our study to the case where the function S to be imbedded is strictly contractive 
at infinity: || ^(00) || < 1. If S is strictly contractive at some finite point iloq, the change of variable 
s — > l/(s — iojq) will make it contractive at infinity and such a transformation preserves rationality 
and the McMillan degree while mapping II + onto itself, hence our results carry over immediately 
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to this case. But if S is strictly contractive at no point of the imaginary axis, then our method of 
proof runs into difficulties and the answer to the minimal degree symmetric inner extension issue 
will remain open. To recap, we pose the following problem: 

Given a p x p symmetric rational Schur function which is strictly contractive at infinity, what is 
the minimal McMillan degree of a 2p x 2p inner extension S of S which is also symmetric : 



S21 



S12 
S 



Su 



Su, S21 



Sf 2 , S(iu)S*(icj) = l2 P , oj € 



(19) 



4.1 Inner extensions of the same McMillan degree. 

Our point of departure will be the solution to the Darlington synthesis problem for a rational 
function in terms of realizations. Let 5 be a Schur p x p function which is strictly contractive at 
infinity and let 



S 



A 


B 


C 


D 



(20) 



be a minimal realization of S of degree n. The strict contractivity at infinity means that I p — D*D 
and I p — DD* are positive definite. Therefore we may set 



A = A + BD*(I p - DD^C, 
B = B{Ip- D*Dy 1/2 , 
C = (I p - DD*)-^ 2 C, 



and subsequently we define: 



H 



-A* 
BB* 



-C*C 
A 



(21) 
(22) 
(23) 

(24) 



Lemma 1 Assuming S is a Schur function strictly contractive at infinity given by \20\). then the 
matrix H defined in \24\) is a dynamics matrix of (I p — 55*) -1 . Furthermore H is Hamiltonian, 
i.e. 



H* 



Ip ' 




Ip ' 


-Ip 




-Ip 



H. 



Proof. By definition 



then, from ([13 



S* 



-A* 


-C* ' 


B* 


D* 





( -A* 





C* 


I p - SS* = 


BB* 


A 


-BD* 




\ DB* 


C 


Ip - DD* 



and if S is strictly contractive at infinity the inverse of I p — DD* is well defined. Then from ([T 
we have 





' -A* - C*A t DB* 


-C*A t C 


C*A l 


(i p - ss*)- 1 = 1 


BA r B* 


A + BD*AiC 


-BD*A t 




-A t DB* 




Ai 
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where = (I p — DD*)~ l and A r = (I p — D*D)~ 1 , whose dynamics matrix is none but H. 
Finally, it is easy to check from the definitions (12 ID - ([23]) that H is a Hamiltonian matrix, i.e. that 
the partition of H defined in ([24j) satisfies H{ 2 = H12, = H21, an d H 22 = — H 



ii- 



□ 



Remark 4.1 The Hamiltonian character of H implies that it is similar to —H*. In particular 
the eigenvalues of H are symmetric with respect to the imaginary axis, counting multiplicities. It 
must also be stressed that the realization of (I p — SS*)~ l given in the proof of Lemma{l\ may not 
be minimal. Because the McMillan degree is invariant upon taking the inverse, the realization in 
question will in fact be minimal if, and only if, the McMillan degree of SS* is the sum of the 
McMillan degrees of S and S* . This will hold in particular when no zero of S is a pole of S* 
[BGl^, in other words if no zero of S is reflected from one of its poles. Hence the characteristic 
polynomial of H plays in the matrix-valued case the role of the polynomial \x given by U6\) in the 
scalar case (compare condition (ii) after (fi 6\)). 



The (not necessarily symmetric) inner extensions of S that preserve the McMillan degree are char- 
acterized by the following theorem borrowed from [GR| . Actually, theorem 4.1 in |GR| describes 
all the rational unitary (on the real line) extensions of a (non necessarily square) rational matrix 
function which is contractive on the real line and strictly contractive at infinity. The next theorem 
essentially rephrases this result in our right half plane setting dealing with square and stable matrix 
functions. 

then all 
(25) 



Theorem 1 If S given by \20\) is a Schur function which is strictly contractive at infinity, 
(2p) x (2p) inner extensions 

Sn £12 

of the same McMillan degree as S are given by 

' Ut 



U 2 








I n 



where U\ and XJ% are arbitrary unitary matrices and where Sp is given by 





( A 


Bi 


B 


s P = 




D n 






V c 


D21 


D 



with 



D 



21 



a/2 



-D*, 



(Ip-DD*) 1 / 2 , D 1 2 = (I P -D*D) 1 

Ci = -(I p - D*D)- 1/2 (B*P- 1 + D*C), 
Bi = -(PC* + BD*)(I P - DD*)- 1/2 , 
and P is a Hermitian solution to the algebraic Riccati equation: 

TZ(P) = PC*CP + AP + PA* + BB* = 0. 

The map P — > Sp is a one-to-one correspondence between the Hermitian solutions to ([57 
inner extensions of degree n of S whose value at infinity is T> defined in HZ 



(26) 

(27) 

(28) 
(29) 
(30) 

(31) 
and the 



Remark 4.2 Note that [GR, thm. 3.4-] guarantees, under the assumptions of Theorem{J\ that all 
Hermitian solutions of \°31\) are invertible and positive definite since S(s) is stable. 
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4.2 Relation to spectral factors. 



We say that a stable p x p matrix-valued function Sl (resp. Sr) is a left (resp. right) spectral 
factor of I p - SS* (resp. I p - S*S) if S L S* L + SS* = I p (resp. S* R S R + S*S = I p ); such a factor 
is called minimal if it is rational and if the block rational matrix (Sl S) (resp. (S^ S T ) T ), whose 
McMillan degree is at least the degree of S, actually has the same McMillan degree as S. This is 
equivalent to require that (Sl S) (resp. (S^ S T ) T ) has a minimal realization whose output (resp. 
input) and dynamics matrices are those of a minimal realization of S. In particular Theorem Q] 
implies that, for any inner extension of S having the same McMillan degree, S21 (resp. Si 2) is a 
minimal left (resp. right) spectral factor of I p — SS* (resp. I p — S*S). 

The corollary below is essentially a rephrasing of the theorem in terms of minimal spectral factors, 
compare [Ml|F3l|W]. Observe that substituting (pgj) - (|50]) in fl27J) yields 







' D u 


D 12 ' 




' B* ' 


c 




. D 2 i 


D 




B* 



p 



-1 



while a similar substitution in (j31j) yields 

AP + PA* + B X B\ + BB* 



0. 



(32) 



Note that (j32|) has a unique solution since A has no purely imaginary eigenvalue, and that this 
solution is necessarily Hermitian positive definite by the controllability of [A B]. 



' A 


B ' 


C 


D 



be a minimal realization of a Schur function strictly contractive 



Corollary 1 Let S 

at infinity, and define Du, D21, and D\2 as in (28\) . To each minimal left spectral factor S21 of 
I p — SS* with value D21 at infinity, there is a unique inner extension of S of the same McMillan 
degree, whose lower left block is S21 and , with value at infinity: 

D n D l2 
D21 D 



V 



If we put 



[S21 S] 



A 


Bi 


B 


C 


D21 


D 



then this extension is none but Sp given by \2_ 
tively, one also has 



(33) 



where P is the unique solution to \3£fy . Alterna- 





r a 


Bi 


B 


S P = 


- (D U B* + D 12 B*)P- 1 


Du 


D 12 




c 


D21 


D 



(34) 



Proof. Suppose we have an inner extension S of S, with the stated properties, given by (I25p . 
Then, from (f3"2"T) - (f3"3|) . the block S21 uniquely defines P thus also S by Theorem [TJ Conversely, 
let S21 be a minimal left spectral factor of I p — SS* with value D21 at infinity. Then we have a 
realization of the form ([33]) for [S21 S] and we may define P through ([32]) . Using (fT3|) . we obtain 
for I p — [S\2 S][Si2 S]* the following realization 



Ip — [S21 S][S2i s] 



-A* 





C* 


BiB{ + BB* 


A 


-B x D* 2l - BD* 


D21B* + DB* 


C 


Ip - D 21 D* 21 - DD* 
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Performing the change of basis defined by 
realization to be 

I P -[S 2 i S][S 21 S}*-- 



Ip 





- p 


h 



using fl28[l and (|32|) . we find another 



-A* 





c* 







-PC* - B1D21 - BD* 


CP + D 2 iBl + DB* 


c 






But the rational function under consideration is identically zero by definition of S21, hence by the 
observability of [C A] we get in particular —PC* — B\D\\ — BD* = which yields 



C = -[D 2 i B]*p-\ 



(35) 



Now, put Ci = -[D u D 12 ][Bt Bfp- 1 and let S be defined as in (|25p . Starting from the realization 



of I2 P — SS* provided by (I13p and performing the change of basis defined by 



P 







In 



computation entirely similar to the previous one shows that this is the zero transfer function, that 
is, S is inner. Moreover, (|35h shows it is an extension of S whose lower left block is S21, and clearly 
it has the same McMillan degree and value T> at infinity. Finally, it is straightforward to check that 
this extension satisfies (|34H . □ 
An inner extension Sp of S with the same McMillan degree and value T> at infinity is thus completely 
determined by the choice of a minimal left spectral factor of I p — SS* . Of course a dual result holds 
true on the right, namely the extension is also uniquely determined by a right minimal spectral 
factor S12 of I p — S*S. In what follows, we only deal with inner extensions having value T> at 
infinity, which the normalization induced by (|26p on letting U\ = U 2 = I p there. 
Let P be a solution of the Riccati equation (|31j) . Then the matrix H defined in (|24h satisfies the 
similarity relation 



Ip 





H 









-p 


h - 


P 


h _ 





-(A + PC*C)* -C*C 

A + pc*c 



Thus 



Xh(s) = xz(s)x-z*(s), 
where xm denotes the characteristic polynomials of M and where we have set 

Z = A + PC*C. 

Moreover, since S is strictly contractive at infinity, S21 is invertible and, in view of (I14p . S. 
the dynamics matrix 



21 



(36) 

(37) 

(38) 
has 



A - BxD^C = A + BD*(I P - DD*)~ 1 C + PC*(I P - DD*y 1 C = A + PC*C = Z. (39) 

Likewise, S^ has the dynamics matrix 

A - BD^d = -P{A + PC*C)*P- 1 = -PZ*P- 1 . (40) 

This way the extension process is seen to divide out the eigenvalues of H between the inverses of 
the left and right spectral factors of I p — SS* and I p — S*S respectively. 
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It is a classical fact that there exists a natural partial ordering on the set of Hermitian solutions to 
(|3ip . namely Pi < P2 if and only if the difference P2 — Pi is positive semidefinite. It is well-known 
(see [LRU sect. 2. 5.] or [LR2) ) that there exists a maximal solution P and a minimal solution P of 
(p?Tj) : P is the unique solution for which a(A + PC*C) C LT + , while P is the unique solution for 
which a(A + PC*C) C II , where a(M) denotes the spectrum of M. The left spectral factor §21 
associated with P is called the outer spectral factor. Its inverse is analytic in II + . 

Proposition 1 Let S be a Schur function which is strictly contractive at infinity. Let P and P be 
two distinct solutions of LSI]) and 





S12 


1 sp — 


Su 


S12 


S21 


s 




S21 


s 



the inner extensions of S associated to them via Theorem d 

Then, the matrix Q = S^i 1S21 is well-defined and unitary on the imaginary axis. Its McMillan 
degree coincides with the rank of P — P and Q is inner if and only if P > P. 

Proof. Assuming that S is strictly contractive at infinity, S^i is invertible and we may define 
Q = S21 £21. We have that 



{ \{ \ * n-l n q* q— * q— 1(T nn*\n- * q — 1/C C* \ C-* T 

W — D 2 i 021»-'21' : '21 — '-'21 l J P — 01:1 ) D 21 ~ D 21 W21'J21/' 3 21 — J P' 



(41) 



so that Q is unitary. A realization of Q can be computed from the realizations of S21 and S21 of 
Theorem [H using flT3]) and CO}: 



Q 



A - B 1 D^C 




(c 


Bi 


-D21C 


D2I ) 




D21 



A - B x D 2 lC 


B x D 2 lC 


B x 





A 


B! 


-D£C 







Applying a change of variables using T 



we get 





(A- BxD^C 





B 1 -B 1 


Q = 




A 


Si 




V -D£C 








A - B x D 2 lC 


B 1 -B l 


-D21C 





(42) 



From ((30]) we draw B x - B x = (P - P)C*D 2 ^ and 

Q 



z 


{P-P)C*D^ 


-D21C 


Ip 



(43) 



where Z is given by f)38[) . 

Set r = P — P. Since V is Hermitian, the singular value decomposition can be written 



V 





r 



v* 



(44) 
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where V is unitary and Tq is real and diagonal. Let d x d be the size of Tq, and partition V as 
V = Wi V 2 ] where V\ is of size n x (n — d) and V2 of size n x d. Note that the columns of Vi span 
the kernel of T. In another connection, it holds that 

K(P) - K{P) = AT + ri* + PC*CP - PC*CP 

= ZT - PC*CT + TZ* -TC*CP + TC*CP + PC*CT 
= ZT + TZ* + TC*CT, 

which is zero since both P and P are solutions to the Riccati equation (|3ip . Thus T is a solution 
to the Riccati equation 

ZT + TZ* + TC*CT = 0. 



Partitioning V* ZV into 



we can rewrite (1451) as 



v*zv 



Z\l Z\i 

Z21 Z22 



(45) 
(46) 





Z\2 




" 


+ 


" 




" z n 


Z\2 


+ 


" 


v*c*cv 


" 


Z21 


Z22 




r 


r 




Z21 


Z22 


r 


r 



0. 



that is 

Therefore 
and 







Z12T0 



TqZ1 2 Z22T0 + TqZ 22 + TqV 2 *C*CV2Tq 
Z 12 = 



Z22T0 + TqZ 22 + TqV 2 C CV2T0 — 0. 
Using V as change of coordinates in the state space, we obtain a new realization for Q: 



Q 



v*zv 


D 2 lV*TC* ' 


CVD 21 ' 





Z U 








Z21 


Z22 




D^CVt 


D 21 CV 2 





(47) 
(48) 

(49) 



where we used (|47p . This readily reduces to 

Z22 



Q 



DZrCV 2 



T V 2 *C*D 2 l 



(50) 



and by (|48p we also have (recall Tq is invertible) that r o 1 solves the Lyapunov equation 

z^o 1 + 1 z 22 + v 2 *c*cv 2 = 0. 



(51) 

From (|50p . it is clear that the degree of Q is at most the size of Tq which is the rank of T = P — P. 
We claim that the realization (|50|) is minimal. Indeed, it is observable because we started from 
the observable realization ([4*2]) and we just restricted ourselves in step ([4*9]) - ([50]) to some invariant 
subspace of the dynamics matrix in the state space. To check reachability, we use Hautus's test 
that no nonzero left eigenvector x T of Z22 associated, say, to some eigenvalue A, can lie in the left 
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kernel of TqV^^D^i", for then (|48p and (|23p together imply that x t Tq, which is nonzero as To has 
full rank, is a left eigenvector of Z% 2 associated to —A. Thus FqX would be an eigenvector of Z22 
and by construction it lies in the kernel of CV2, contradicting the observability of (|50p . This 
proves the claim, to the effect that the degree of Q is in fact equal to the rank of T = P — P. Now, 
from classical properties of solutions to Lyapunov equations |BGRl th.6.5.2], the number of poles 
of Q in II + is equal to the number of negative eigenvalues of Tq 1 solving (15ip . □ 

Remark 4.3 A similar result holds true for the right inner factors: the matrix R = S12S12 is 
unitary on the imaginary axis and inner if and only if P < P. 

The outer spectral factor play an important role in what follows, due to the fact that any (not 
necessarily minimal) left spectral factor, say, a of I p — SS* can be factored as a = 021Q where Q 
is inner. Indeed, the strict contractivity of S at infinity entails that a is invertible as a rational 
matrix, and then computation (j41|) with S21, S21 replaced by S21, o shows that the rational matrix 
Q = S 2 i o~ is unitary on the imaginary axis. In particular it cannot have a pole there, and since it 
is analytic in II + as a is stable and S21 outer, we conclude that Q is stable thus inner, as desired. 



4.3 Inner extensions of higher degree 

We shall be interested in inner rational extensions where we allow for an increase in the McMillan 
degree, and we will base our analysis on the next proposition. For the proof, we need a notion 
of coprimeness: two inner- functions L\, L2 are right coprime if one cannot write L\ = G\J and 
L>2 = G2J with G%, G 2 , J some inner functions and J non-constant. It is well known [Fl} [BOl 
that this is equivalent to require the existence of two stable rational matrix functions X\ and X2 
such that the following Bezout equation holds: 

X 1 (s)L 1 (s)+X2(s)L 2 (s) = I p . (52) 

Left coprimeness is defined in a symmetric way. 



Proposition 2 All rational inner extensions of a Schur function S, contractive at 00, can be 
written on the form 



' L 





s P 


' R 
















(53) 



where L, R and Sp are inner, and Sp is an extension of S at the same McMillan degree, obtained 
from a solution P of the Riccati equation $31}) . 



Proof. Let 



<Sll S\2 

S21 s 



be the inner extension of S of degree n associated with the minimal solution P of f)31 [) . so that S21 
is the outer left spectral factor of I„ — SS* . For an arbitrary rational inner extension 



o~u a 12 

0~21 S 
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of S, we mechanically obtain since EE* = SS* = Ii v that 



S21S21 = Ip — SS* = o"2iO"2i and 



Ol2iJ12 = lp — O O = (T 12 C r 12- 



Therefore, computing as in (jlTT) . we can write a 21 = S21R and 012 = LSvz where R = S^o^i 
and L = cx^S^ 1 are rational matrices that are unitary on the imaginary axis. By the discussion 
after Remark 14.31 R is inner. Next, using again that EE* = I2 P , we get an = — a^S*^* = 
— LS12S* S21 R = LS\\R, and therefore 



" L 







" £11 


£12 




' R 








I P . 




£21 


S 










If L is inner, we have finished because fj53|) holds with P = P. Otherwise, being a unitary rational 
function, L can be factored as L = L\L\ where L\ and L2 are right coprime inner functions, 
so that (|52p holds for some stable transfer functions X%, X2. In fact, the existence of such a 
factorization follows from the so-called Douglas-Shapiro-Shields factorization [DSS] as carried over 
to matrix- valued strictly non cyclic functions in [Flj . see e.g. [BOl] for a detailed discussion of the 
rational case in discrete time that translates immediately to continuous time by linear fractional 
transformation. 

From (j52j) we deduce that L^s) = Xi(s)L(s) + ^(s), therefore 



(54) 



T * 









S12 




' R 







' Xi 





s + 


" x 2 





s 


' R 








Ip _ 




5*21 


S 







I P _ 







I P _ 







Ip _ 





I P _ 



is stable. In particular, if we put S 
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^2^12) we see that 



S12 




T * 







S\2 


s 







Ip _ 




S 



(55) 



is stable. Now, multiplication by a rational function analytic in II (including at infinity) followed 
by the projection onto stable rational functions (obtained by partial fraction extension) cannot 
increase the McMillan degree; this follows at once from Fuhrmann's realization theory |F1} IBOT] , 
Therefore, we deduce from ([55]) that the degree of [Sf 2 S T ] T is at most the degree of [S[ 2 S T ] T . 
But the latter is equal to n for Corollary Q] as applied to S T implies that S\2 is a minimal right 
spectral factor of I p — S*S. Hence the degree of [Si 2 S T ] T is n and Si 2 is again a minimal right 
spectral factor of I p — S*S. Thus by (the transposed version of) Corollary [H there exists an inner 
extension of S of McMillan degree n associated with S12: 



S P 



Sn S12 
S21 s 



As seen after Remark 14.3] one has S21 = S21R2 for some inner matrix i?2- Moreover, by the inner 
character of Sp and Sp and in view of (f55j) . we get 



50 o* c — * T * c c* c — * v> — * r * O D _ 

11 — — <5i2<5 <J2i — — iv2<Ji2<J <J21 2 — -^-'2 1 1 -^2 



Altogether, this implies 



Sn 


S12 











T * 
n 2 







Su 


S12 


S21 


s 







Ip . 







Ip . 




S21 


s 



(56) 
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Next, we contend that the inner functions Sp and diag[i?2)^p] are right coprime. Indeed, if R c is 
a right common inner factor, we get in particular R c = diag [R2, Ip]J* for some inner J, so the last 
p rows of both J and J* are analytic in n + , i.e. they are constant. This entails 



Rr = U 



R3 

I n 



(57) 



for some constant unitary matrix U. But then SpR*U is inner since R c divides Sp, and it is an 
extension of S by (15T|) . Taking determinants, we see that this extension has McMillan degree equal 
to deg5p — degi?3. As this degree is at least equal to that of S which is also that of Sp, we conclude 
that i?3 thus also R c are constant, which proves our contention. 
By the coprimeness above, there exist stable rational matrices X and Y such that 



X 



Sn 


S12 


+ Y 


' R2 





S21 


S 





Ip - 



»2p- 



(58) 



From ([5l|) - (f5"6"]) . we see that the product Sp diagfit^i?, I p ] is stable. Therefore, right multiplying 
(|58p by diagfit^i?, Ip], we deduce that the latter is also stable hence R\ := R^R is inner. Finally, 



" Li 


s P 


' Ri 





Ip. 


.° i P . 






is indeed of the form (|53|). as wanted. 



□ 



5 Symmetric unitary extensions 

We assume from now on that the Schur function S is symmetric i.e. S = S T and we consider a 
symmetric realization, 



S 



A 


B 


C 


D 



, A = A T , B = C T , D = D T (59) 



Such a realization always exists thanks to Theorem 5 in [FHJ. 

It follows that the matrices A,BB*,C*C defined by ([21]), (J25D and ([23]) satisfy 

A = A T 
BB* = (C*C) T . 



(60) 



If the extension Sp associated to some solution P of f|31 j) raa Theorem [T] is symmetric, then the 
matrix Z defined in (|38|) must be similar to — Z* since they are, by the computations (|39p and 
([I0|) . dynamics matrices of S^i an d respectively. Therefore, in view of ([37]) . the characteristic 
polynomial \H must be of the form 7r(s) 2 . As this may not be the case, a symmetric inner extension 
of S preserving the McMillan degree may well fail to exist. However, as we shall see {cf. also [AVJ), 
symmetric inner extensions of higher degree do exist. We will first give a simple method to construct 
(possibly unstable) symmetric unitary extensions of S using the inner extensions Sp provided to 
us by Theorem Q] 

As the Riccati equation f|31|) does admit a Hermitian solution and since the pair (A, C) in (|2ip - 
(|23p is observable, as follows immediately from Hautus's test on using the observability of (A, C), 
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the partial multiplicities (i.e. the sizes of the Jordan blocks) of the pure imaginary eigenvalues 
of H (if any) are all even, see [LRU th.2.6] or [LR2} th. 7.3.1]. Let 2no be the dimension of the 
spectral subspace of H corresponding to all its pure imaginary eigenvalues. Using (|36p and grouping 
together the roots of even multiplicity, we can write the characteristic polynomial of H in the form 



Xh(s) = vr(s) 2 x+(s)x K 



(61) 



where vr(s), Xk ( s )> an d Xk ( s ) are polynomials in s, and where xt has k simple roots in II + while 
X~ = (xt)* ■ Then k is the number of distinct eigenvalues of H in II + with odd multiplicity, and 
7r(s) has degree greater than or equal to uq by what precedes. 

Proposition 3 Assume S is a symmetric Schur function strictly contractive at infinity, and let 
l5Pj) be a symmetric realization. For P a Hermitian solution of [31\). let further 



S P 



S21 



S12 

s 



be the inner extension of the same degree associated by Theorem [IJ 



Then, the matrix Q = S 21 Sf 2 * s well-defined, unitary on the imaginary axis, and 



Sn 


S12 




' Q 







Dil0 2 l °12 


Sl2 


S21 


S 











oT 
°12 


S 



is a symmetric extension of S which is unitary on the imaginary axis. The McMillan degree of Q 
coincides with the rank of P~ T — P , and Q is inner if and only if P~ T — P is positive semi-definite. 
In this case, Sp is inner and its McMillan degree is degS + degQ. 

Moreover, deg Q is in any case greater than or equal to n, the number of distinct eigenvalues of H 
in II + with odd multiplicity. 



0, and that S p -t 



dp. 

^21 ^12 ^ S 



Proof. From §0$) we see that TZ{P) = if and only if TZ(P~ T ) 
Therefore, we may appeal to Proposition Q] with P = P~ T and S21 = Sf 2 . Thus, Q 
unitary on the imaginary axis, its degree d is equal to the rank of P~ T — P and Q is inner if and 
only if P~ T > P. The computations in the proof of Proposition Q] apply, so if (|44p is the singular 
value decomposition of T = P~ T — P and Z is defined by (I38p . we get from (I46p and (I47p that 



v*zv 



v x zv x v x zv 2 

V 2 *ZVx V 2 *ZVx 












Z21 


^22 



(62) 



On the other hand, if we set Z := -PZ*P~ 1 , it follows from flMj and (|38D that 



Z 



-pz*p- 



-p(a* + c^pyp- 1 = (ap + bb*)*p- 



But since the columns of V\ span the kernel of T, we get V£(P — P T ) - 
P T Vi = P'Wi, hence in view of (JB3J) and dSDJ 

ZVi = (i + BB*P T )V 1 = Z T Vi. 

As V is unitary, we have thus arrived at a similarity relation of the form: 



A + BB*p-\ (63) 
= that transposes into 



V T ZV 



V^ZVi V?ZV 2 
V 2 T ZVi V 2 T ZV 2 







Z\2 
Z22 



(64) 
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Now, since Z is similar to — Z*, we see from ([62]) . ([64]) . and ([37]) that x| n ( s ) must divide 
Consequently the number of rows (or columns) of Z\\, which is the dimension of the kernel of 
P~ T — P, cannot exceed the degree of 7r(s) in f)61 [) . namely n — k. Therefore rank(P _T — P) > k, 
as announced. □ 

Remark 5.1 Note that degQ = 0, i.e. Q is constant, if and only if P~ T = P. 



6 Minimal symmetric inner extensions 



Lemma 2 Let P and P be the minimal and the maximal solution to the Riccati equation \31)) 
associated with a symmetric Schur function strictly contractive at infinity. Then 

p-T _ p 

Proof. Since P is the minimal solution, for each solution P we have that P > P. By symmetry, 
P is a solution if and only if P~ T is a solution and moreover P~ T > P~ T , so that P~ T must be 
the maximal solution. □ 
A symmetric inner extension of the symmetric Schur function S may now be obtained as follows. 



Proposition 4 Let 



Su 
S21 



S12 
S 



be the extension associated with the minimal solution P to the Riccati equation (El 
extension T,p given by Proposition^ 



The symmetric 





S12 




' Q 





S21 


s 










with Q = S, 



-1 qT 

21 °12> 



is inner and has degree 2n — tlq. 



Proof. Since P is the minimal solution of the Riccati equation, we have that P T — P > 0, hence 
by Proposition [3] the extension Sp is inner and has degree n + d where d is the rank of P~ T — P. 
By Lemma[2]P _T = P, and if we let ker(P — P) = J\f it is a classical fact (see e.g. [LRTJ th.2.12] or 
[LR2J) that J\f is, for any solution P to (f3Tjh the spectral subspace of Z = A + PC*C corresponding 
to all of its pure imaginary eigenvalues, and that it has dimension uq. Therefore deg Sp = 2n — tlq, 
as desired. □ 
Finally, we shall construct from Ep a symmetric inner extension of S of minimal degree. For this, 
we first establish a factorization property which strengthens the Potapov factorization (J7J) in the 
symmetric case. 

For Re £ > 0, recall from ([5]) and © the definitions of b^ and B^. Let U be a unitary matrix whose 
first column is a unit vector mGP, then 



B^ u = UB^s)U* = L P + (b^s) - l)uu* 
is an elementary Blaschke factor with the following properties (compare |D4|, chap.l]): 



detP ?jU (s) 



l)uu* 



(65) 



(66) 
(67) 
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Note that any elementary Blaschke factor can be written in the form B^ )U V for some unit vector u 
and some unitary matrix V . 

Lemma 3 For F a p x p symmetric complex matrix, there exists a unitary matrix U such that 

F = UMJ T (68) 

with A a non-negative diagonal matrix. If F has rank k, the p—k first columns ofU may be construed 
to be any orthonormal basis of the kernel of F. 

Proof. The existence of (I68p is proved in [HJ} cor. 4. 4. 4] under the name of Takagi's factorization. 
Permuting the columns of U, we may arrange the diagonal entries of A in non-decreasing order and 
then the first p—k columns form an orthogonal basis of KerF. Right multiplying U by a block 
diagonal unitary matrix of the form diag [V, Ik] we may clearly trade that basis for any other. □ 

Lemma 4 Let T(s) be a symmetric inner function and £ G II + a zero ofT(s). Suppose that there 
exists a unit vector u such that the interpolation conditions 

T(£,)u = (69) 
u T T'(£)u = (70) 

are satisfied. Defining B^ u (s) as in \65}) . then, 

R{s) = B^ u {sr T T{s)B^ u {s)-\ 

is analytic at £ and thus a symmetric inner function of degree N — 2. 

Proof. We give a proof of this result which follows Potapov's approach to the multiplicative 
structure of J-inner functions [PJ. For simplicity, we use Landau's notation 0{s) k for the class of 
(scalar or matrix- valued) functions f(s) such that ||/(s)||/|s| fc is bounded for \s\ sufficiently small. 
We also put ||/(s)|| for the operator norm of f(s). Write the Taylor expansion of T(s) about £ as 

T(s) = T(0 + (s - 0T'(0 + 0(s - 2 , 

and apply Lemma [3] to obtain Takagi's factorization of T(£) in the form 

U T T(0 U = diag(0, . . . , O.pi, . . . , p r ), (71) 

for some unitary matrix U whose first column is u. Next, define a matrix- valued function T\ by 

Ti(s) := T(s)B^ u (s)- 1 

= TftB^s)- 1 + (s- OT'iOBtAsy 1 + 0(s - o 

= (CT) T diag(0, . . . , 0, Pl , . . . , Pr ) B^(s)- l U* + (s - ^(OB^i 8 )' 1 + °( s " 

= T(0 + (s- OT'iOB^uis)- 1 + 0(s - 

where we have used (f65|h Clearly, Ti(s) is analytic about £ and, since B^ ;U (s) has degree 1, T\(s) 
is an inner function of McMillan degree N — 1. Now, using (j66|) and the symmetry of T(£) which 
implies m t T(£) = 0, we get the following interpolation condition for T\(s): 

u T Ti(£) = u T T(£) + 2Re £u T T'(£)uu* = 0. 
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Applying what precedes to Ti(£) T , we see that R(s) := Ti(^) T B^ u (s) 1 is analytic at £ and thus 
an inner function of degree N — 2. Finally, we can write 

R(s) = B^ u (s)- T T(s)B ( , a (s)- 1 , 

which achieves the proof. □ 

Proposition 5 Let T(s) be a symmetric inner function of degree N and suppose that T(s) has a 
zero £ G II + of multiplicity greater than 1. Then there exists an elementary Blaschke factor of the 
form B(s) = B^ u , where u G C p is a unit vector in the kernel o/T(£) ; such that 

R(s) = B(s)~ T T(s)B(s)- 1 

is analytic at £. Thus R(s) is inner, symmetric, and it has McMillan degree N — 2. Moreover, if 
V is a 2-dimensional subspace of the kernel o/T(£), we may impose in addition that u G V. 

Proof. Write the local Smith form of T at £, namely (jlip where S is replaced by T and k is equal 
to p. Assume first that the kernel of T(£) has dimension 1 over C. Since T is analytic at £ and £ 
is a zero of multiplicity at least 2, the partial multiplicities must satisfy Vj = for 1 < j < p and 
Up > 2. Therefore, if e p is the last element of the canonical basis of C p , the vector uq = -F _1 (£)e p 
spans the kernel of T(£) and the C p -valued function = F~ 1 (s)e p is such that T(s) cj>(s) has a 
zero of order at least 2 at £: 

T( S )0( S ) = O(s-£) 2 . (72) 
Let us write the Taylor series of T(s) and about £: 

T( S ) = T(£) + ( S -£)T'(£) + 0( S -£) 2 , 
0(s) = u + ui(s-£) + O(s-£) 2 , 

for some u\ G C p . Then 

T{s)4>{s) = T(Ou + (s- £)(T'(£)n + r(£W) + 0(s - £) 2 , 

and we must have 

T(Ou = 
T'(£)n + T(£)m = 0. 

By symmetry the first equation implies UqT(^) = and then the second one yields UqT'(^)uq = 0. 
Since «o 7^ 0, we obtain the desired result from Lemma H] with u = uq/\\uo\\. 

Assume next that the kernel of T(£) has dimension at least 2 over C and let V be a 2-dimensional 
subspace. Two cases can occur: 

(i) There exists a non-zero v G V such that T{s)v vanishes at £ with order at least 2. On grouping 
the partial multiplicities in such a way that er,- = for 1 < j < jo, o~j = 1 for jo < j < j\, and 
Oj > 2 for ji < j < p, we deduce that in the decomposition 

jo jl v 

F{s)v = ^2<Pj(s)ej + <A?( s ) e i + Y <M s ) e i 

i=i i=io+i j=ji+i 



19 



the functions <pj vanish at £ with order at least 2 for 1 < j < j\ and at least 1 for jo < j < Ji- In 
particular F'(^)v lies the span of the for jo < j < p which is the kernel of T(£)i ?_1 (£). Therefore 
if we set cfi(s) := F~ 1 (s)F(S i )v, we get an analytic function about £ such that T(£)(f)(!;)v = and 

(T0)'(O = T'(Ov - T&F-^OF'iOv = + = 

so that (|72p holds. The result now follows as in the previous part of the proof with uo = v. 

(ii) Each non-zero v £ V is such that T(s)v vanishes at £ with order 1. By Takagi's factorization 

t/ T T(£)[/ = diag(0,...0,p 1 ,p 2 ,...p r ), 

for some unitary matrix U, where < r < p — 2 and < p\ < p 2 < • • • p r - We may assume that 
the first two columns of U form a basis of V. Define a matrix-valued analytic function about f; by 

Q(s) := [/ T r( S )£7 = diag(0, . . . , 0, p x , p 2 , . . . , p r ) + (s - 0Q'(0 + 0(s - £) 2 - 

The matrix is symmetric, and its 2 x 2 left upper-block is invertible. Otherwise indeed, 

there would exist a non-zero y = [2/1,2/2, 0, ... , 0] T £ C p such that T(s)Uy vanishes at £ with order 
2, but this contradicts our standing assumption since Uy £ V by the choice of U. 
We now exhibit a non-zero vector w which satisfies the interpolation conditions 





0. 



(73) 
(74) 



To meet ([73jl . we put = (x, 0, . . . , 0) T with x G C 2 so that w T Q'(^)w = x T R\(^)x. By LemmaH 
we may write Uf Ri(£)Ui = A, with f/i is a 2 x 2 unitary matrix, where A = diag[A 2 , A|] is strictly 



positive. Define now v 



[ A2 i\i ] . Obviously 



v T Av 



1 



Af + A^ 



[A 2 , -iAi] 



hence setting x = U\v we have th 
let u = U[w T , 0] T , we get u G V and putting £?(s) 



" A? " 




A 2 


. A 2 _ 




— z'Ai 


x T i?i(0^ 




v T Av = 



B*(s) T(s)B*(s) is analytic about £, which achieves the proof. 



0. 



0, as desired. Finally, if we 
l)uu* we see from Lemma 0] that 

□ 



Theorem 2 Let S be a p x p symmetric Schur function of McMillan degree n which is strictly 
contractive at infinity. Let further H be the state characteristic matrix defined by \21$ for some 
minimal symmetric realisation (A, B,C, D) of S. Assume that xh has exactly k distinct roots in 
II + with odd algebraic multiplicity. Then, S has a 2p x 2p symmetric inner extension of McMillan 
degree n + k, and this extension has minimal degree among all such extensions of S. 

Proof. We apply Proposition [5] to the symmetric inner extension of Proposition [3J 



Q 
/„ 



in which P is the minimal solution to (|3ip . Recall that P is the only solution to the Riccati equation 
such that a(A + PC*C) C S~ (see section E£2} , or equivalently a(A* + C*CP) C n+. By fl36]), 
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(|39p . and Proposition U the eigenvalues of H in II - , which are (counting multiplicity) n — no in 
number by (|6ip . are precisely the poles of Q. Likewise the eigenvalues of H in II + are the zeros of 
Q, that are reflected from its poles across the imaginary axis with corresponding multiplicities. 
Let £1 £ n + be one of these with multiplicity strictly greater than one. By proposition \5\ there 



exists an inner function B = B^ lU such that T\ := B UpB is inner of degree 2n — uq — 2. 
Clearly T\ is symmetric and moreover, since 



Ker 



C Ker£p(£i), 



it follows from the proposition that u may be chosen of the form [v? , T ] T (each block being of size 
p) in which case B is of the form 

" B 
I P 

so that T\ is an extension of S. The matrix H has n — uq = k + 21 eigenvalues in n + , counting 
multiplicities. Thus, we can perform i iterations to obtain an extension of degree 

2n — no — 21 = n + k. 

We now prove that this extension has minimal McMillan degree. Let £ be any symmetric extension 
of S. By Proposition O it can be written in the form 



" L 





s P 


' R 


















with Sp 



S21 S22 



where L, R, are inner, and Sp is an inner extension of S of the same degree. The extension X being 
symmetric, we must have 

(LS12) = S21R S21S12 = RL. 
By Proposition El the degree of the unitary matrix S^S^ cannot be less than k. This yields 

-IcT 



so that 



n < deg S 21 L S{2 = deg RL < deg R + deg L, 
n + k < n + deg R + deg L = deg S. 



□ 



Corollary 2 Let S be a symmetric Schur function of size p x p which is strictly contractive at 
infinity. Then, the following propositions are equivalent 

(i) S has a symmetric inner extension of size (2p) x (2p) of the same McMillan degree, 

(ii) there is a Hermitian solution P to the algebraic Riccati equation i31\) that satisfies PP T = I p , 

(Hi) the characteristic polynomial of the matrix H given by A21\) -(24\) can be written 

Xh(s) = tt(s) 2 (75) 

for some polynomial tt. 
Moreover, all extensions (i) are parameterized by i26)) . where U\ = and P satisfies (ii). 
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7 Completion of a real Schur function 



We now assume that the Schur function S(s) is real that is to say 



S(s) = S(s). 

It thus admit a realization (|20p in which A, B, C and D are real matrices. We have that 



Theorem 3 Let P be a Hermitian solution to h31\). Then, the associated minimal unitary extension 
Sp of S is real if and only if P is real. 

Proof. Formula (|27p for Sp shows that if P is real, then Sp is real. Conversely, assume that Sp 
is real. Formula (|27p shows that 

C(sl n - A)~ l PC* 

is real for real s. This implies that CA k PC* is real for k = 0, 1, ... , and from the observability of 
(C, A) it follows that PC* is real. The equation can be rewritten in the form 

PC*(I P - DD*)~ 1 CP + (A + BD*{I P - DD*y l C)P 

+P(A* + C*(I P - DD*Y 1 DB*) + B(I P - ZTD)-^* = 0, 

and we see now that AP + PA* is real. Write P = P\ + iP2 where Pi and P2 are real matrices. 
Since CP is real, we have CP2 = and P%A* + AP2 = 0. Premultiply the later equation by C to 
obtain CAP2 = 0, and by induction CA k P2 = 0, k = 0, 1, . . . ,. Now, by the observability of (C, ^4), 
P2 = and P is real. □ 
A real minimal realization of a symmetric function may fail to be symmetric. However there exist 
a minimal real realization which is signature symmetric (see [YT] or [F2} th.6.1]), that is to say 

A T = JAJ, B T = CJ, C T = JB, D T = D, 



for some signature matrix 

J 



"In—r 



It follows that A T = JAJ, (BB*) T = JC*CJ, (C*C) T = JBB*J, which implies that if P is 
a solution of (f3T|) . then P = JP~ T J is also a solution and we have Sp = Sp. If Sp is a real 
inner extension of S computed as in Theorem [TJ then an analog of Proposition [3l in which P~ T 
must be replaced by P, allows to construct a real symmetric unitary extension of S with the same 
properties. However, the situation is more involved than in the complex case. Even if the algebraic 
multiplicity of the eigenvalues of H are all even, a symmetric real extension at the same degree 
may not exist. We conclude by illustrating this with an example. 



7.1 An example 

Consider the rational matrix 

S(s) 



/W 

m 
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where /(s) is a strictly proper scalar Schur function of McMillan degree 1. Let 

f(s) = d + c(s — a)~ x c, 

be a symmetric realization of f(s). The matrix- valued function S(s) is Schur and contractive at 
infinity and has the minimal realization 

5(a) = D + C{sl 2 - A)- l C T , with D = dl 2 , C = cl 2 , A = al 2 . 



In this example, all the eigenvalues of H defined in (|24p have even algebraic multiplicity. We 
first construct a complex symmetric extension at the same degree, degree 2. The Riccati equation 
associated with S(s) is (see th.3.1) 

C*CP 2 + (A + A*)P + BB* 

where 



0, 



(76) 



A 



a(l-\d\ 2 ) + c 2 d 



l-|dP 

The equation (|76p can be rewritten as 



/ 2 , 



1 — Ml 



r/2, <?*C 



r/ 2 - 



p2 | 2 Re(A)(l^ 



\d\< 



-P + h 



0. 



(77) 



A Hermitian solution to (|77|) can be diagonalized in the form P = £7diag(A, fJ>)U*, where U is an 
unitary matrix. Note that the eigenvalues of P must be solutions to the scalar Riccati equation 

M 21 



-p + l = Q, 



(78) 



and must be positive real, since /(s) is Schur (positive real lemma). By corollary [21 a complex 
symmetric extension is obtained from a solution P to (|77p which satisfies PP T = I2. Such a 
solution is obtained taking 



V2 



1 



and for A and [i the two distinct solutions to ([78 p which satisfy A/x = 1. The solution is 

A + /i i(A — /i) 



2 



-i(A — /i) A + n 



It is easy to check that P is Hermitian and moreover that PP T = h- 

We now come to the real case and assume that a, c and d are real. The Schur function S(s) has a 
real extension if and only if the Riccati equation (|77j) has a real solution P. Such a real solution 
must be of the form 



P = 



A 








o 1 



where O is a real orthogonal matrix, and A and fi are (positive real) solutions to the scalar Riccati 

equation (fTHj) . Note that 5(s) is Schur if and only if ([75]) has a positive real solution, which happens 

2 

if and only if a < —^d- A symmetric real extension of S(s) at the same degree is obtained from a 
real solution P to (|77p which in addition satisfies PP T = I 2: but this may only happen if P = 1%, 
that is when a = — c 2 /(l — d), in which case the two solutions of (|78p are equal to 1. For example, 
the function /(s) = l/(s + £) is Schur if and only if £ > 1. It has a real symmetric extension of 
degree 2 if and only if £ = 1. 
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